Superconducting rf-wave junctions: The disappearance of the odd ac components 
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We study voltage-biased superconducting planar d-wave 
junctions for arbitrary transmission and arbitrary orientation 
of the order parameters of the superconductors. For a certain 
orientation of the superconductors the odd ac components 
disappear, resulting in a doubling of the Josephson frequency. 
We study the sensitivity of this disappearance to orientation 
and compare with experiments on grain boundary junctions. 
We also discuss the possibility of a current flow parallel to the 
junction. 

PACS numbers: 74.50.+r, 74.20.-z 



tion of the superconducting order parameter with respect 
to the junction interface as explained in Fig. [l| The cases 
a = and a = 7r/4 correspond to d x 2_ y 2 and d xy sym- 
metry, respectively. For the s-wave case, A(0) = A s . 

For conventional superconducting junctions with s- 
wave electrodes the frequency of the non-stationary 
Josephson current is luj = 2eV/h. This means that the 
Josephson current perpendicular to the junction can be 
decomposed into papponents with frequency parts mujj 
(m is an integer) £11113 



J(i) = ^/ m e 



imujjt 



(i) 



Some of the hole-doped high-T c cuprates show fea- 
tures that cannot be explained using an s-wave or- 
der parameter .Era Instead various phase-sensitive exper- 
iments involving Josephson junctions suggest that the 
order parameter changes sign at certain regions of the 
Fermi surface. The outcome of these and other experi- 
ments makes d-wave symmetry a particularly interesting 
candidate for the symmetry of the superconducting order 
parameter in the high-T c cuprates.u 

Hu's prediction of midgap states (MGS) at inter- 
faces/surfaces of d-wave superconductorscl called for a 
reinvestigation of a number of transport problems involv- 
ing d-wave superconductors. This reinvestigation has to 
some extent been carried out, starting. with the normal 
metal-superconductpt. (NS) junctionpj followed Wr4he 
dc Josephson effcctETEl and the ac Josephson effect .Bli3 It 
has been found that the zero-energy MGS influence the 
current-voltage relation (the NS-case and the ac Joseph- 
son effect), the temperature dependence, and the current- 
phase relation (the dc Josephson effect). 

The origin of MGS is due to normal scattering of the 
quasiparticles at interfaces/surfaces. Since the quasipar- 
ticle changes its momentum when scattered, it will in the 
d-wave case experience a different order parameter after 
scattering. When there is a sign change of the order pa- 
rameter after scattering, a bound state is formed at the 
interface/surface of the superconductor B 

In this paper we investigate the ac Josephson effect 
in planar d-wave junctions, calculating the first Fourier 
components of the current. The order parameter A(9) for 
an anisotropic superconductor depends on the angle 9 of 
incidence for the quasiparticle approaching the junction. 
For the d-wave case the order parameter is taken to be 
A(#) = Aq cos[2(6> — a)], where a is the angle of orienta- 



The zero-frequency|-Gpmponent has been studied recently 
for the d-wave case.t2l In this paper we extend the analysis 
to include the nonzero frequency parts of the current. 
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FIG. 1. Layout of the junction. The superconductors can 
be rotated to an angle 01,2- The strength of the barrier can 
be tuned all the way from the ballistic to the insulating case. 



In the s-wave case one finds that all components I m 
in Eq. (|l|) are nonzero although they decrease with in- 
creasing m. This is not necessarily the case when d-wave 
superconductors are introduced in the junction. Consider 
the d-wave junction and let the order parameter of the 
left (right) superconductor be Ai(0) [A2(0)]. For orien- 
tations of the superconductors where Ai(#) = Ai(— 0) 
(the case for s-wave and d x 2_ y 2 orientation of the left su- 
perconductor) and ^2(6) = —A2(—6) (d xy orientation of 
the right superconductor) we find that I m = for odd m 
and therefore the Josephson frequency is doubled. This 
doubling of the Josephson frequency has previously been 
inferred from the stationary Josephson effect Jij where 
one can show that in the situation described above the 
Josephson current is 7r-periodic in the phase rather than 



1 



27r-periodic. This change of period in the phase would 
then change the Josephson frequency accordingly. 

Considering the nonstationary Josephson effect, we 
may calculate the Fourier components and explore the 
sensitivity of the odd ac components to deviation from 
orientations where disappearance of the odd ac compo- 
nents occur. This is the main subject of this paper and 
is not possible to treat within the framework of the sta- 
tionary Josephson effect. 

In passing, we note that the disappearance of the odd 
ac components is not related to the presence of MGS 
itself; even for the ballistic junction (when there are no 
MGS present) the odd ac components disappear in the 
situation described above. 

In our calculation self-consistency of the superconduct- 
ing order parameter is not fulfilled. Recent investigations 
in the tunneling limit have considered self-consistency 
and found that interface states with nonzero energy ap- 
pear along with the zero-energy MGS El In this paper we 
neglect these effects. 

The system that we consider is shown in Fig. [j]. Two 
superconductors are separated by a barrier represented 
by a scattering matrix, which for electrons is 



Se(9) = 



r{6) 



t{9) 



t(9) -r*(6)t{6)/t*{6) 



(2) 



The scattering matrix for holes is Sh{9) = S*(9). The 
normal reflection at the barrier changes the angle from 
9 to 9 — 7r — 9. All barred quantities in this paper are 
related to 9. 

We solve the time-dependent Bogoliubov-de Gennes 
equation for <i-wave superconductors piecewise in each 
region. The details of the method can be found in e. g. 
Ref. |l(| The order parameter is assumed to be 

f Ai(0) , x<-L/2 
A(k,x,t) = < , \x\ < L/2 (3) 

[ A 2 (9)e^ 0+2eVt ^ , x > L/2. 

The overall phase is not important and we therefore 
choose Ai real and let the phase difference between the 
superconductors be c/)q. In Eq. (J3j> it is understood that 
we only consider k- vectors with |k| = hp. 

Solving the scattering problem, we inject an electron- 
like quasiparticle from the left superconductor at energy 
E and angle 9. Since there is a voltage bias over the 
junction, the quasiparticle goes through multiple An- 
dreev reflections (MAR). The barrier is modeled as a 
(5-function potential of strength H with transmission am- 
plitude t(9) = icos9/(icos9 — Z) and reflection am- 
plitude r(9) = Z/(i cos 9 — Z), where Z is defined as 
Z = 2mH/h 2 . Our calculations include normal regions 
to the left and to the right of the barrier, but the length 
L is put to zero at the end of the calculation. Therefore, 
our method-follows closely the one presented by Averin 
and Bardas.Ej However, the same results would have been 
obtained if we-ipd used the method of Bratus, Shumeiko, 
and WendinJliy 



It is convenient to calculate the current in the normal 
region to the left of the barrier. The wave function in 
this region is 



E 



x 2n c 



•/k-x 



(4) 



where the momentum is defined by k = (k x ,k y ) — 
k(cos9,sm9) and k = (—k Xl k y ) = fc(cos sin #). We 
have introduced the energy E n = E + neV , where n is 
an integer. Note that we have separated out the phase 
difference cf>o between the superconductors, making the 
coefficients a, b, c, and d independent of </>o- Matching 
the wave functions of the different regions one gets the 
coefficients a, b, c, and d as written down in Ref. [l^. 

The next step is to calculate the current by inserting 
the wave function into the following formula for the cur- 
rent density 

J ' = ~^ ZZ tmih \2k~Tj Im ^» VU " + V n^ V n} , ( 5 ) 

where (u n ,v n ) is the solution of the BdG equation. The 
sum is over all incoming electron-like quasiparticlcs for 
both negative and positive energies. There are different 
ways of writing the current density. Using the|-cpmplete- 
ness relation for solutions of the BdG equationlla one can 
show that Eq. (||) is equivalcnL-tp-.the current formulas 
written down in previous workHjr 6 H 17 l 

Including also the left-moving quasiparticles we can 
determine the expression for the ac components I m in 
Eq. #: 



I m (V) 



= c m (vV (m0 °- Qm) , m^o 



C m (V) = ^/Al l (V)+Bl(V), a m = arctan 

2 5/2 em l/2£l/2 Ao ^ 



O"0 — Ly 



h 2 



B m (V) 

A m {vy 

(6) 



where D = J d9\t\ 2 cos 9/2 is the transmission probability 
through the barrier averaged over all angles 9, and L y is 
the junction length in the y-directionil3 In Eq. (^J) we 
introduced A m (V) and B m (V) which are the coefficients 
in the cosine and sine expansion of the ac current. The 
coefficients are defined as 



A m{V) = jf, d9cos9 / —tanh 

J -tt/2 7-oo A \2k B T 

x ^ TN T (E)Re {T T (E, 9, m}} , 



Bm{V) = -— / dffcosd / —tanh 

4DJ_ 7r /2 J-oo A \2k B T 

xJ2N T (E)lm{T T (E,9,m)}, (7) 
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where r = + (— ) means right (left) movers, N T (E,9) is 
the bulk superconducting density of states evaluated in 
superconductor 1 (2) for r = + (— ), and 



T T (E,9,m) 



( fl 2n+2m,T a 2n,T d 2n+2mT d 2r , 



n— — oo 



+ ^2n+2m,r^2n,T c 2n+2m,r C 2n,r) • (8) 



The phase m<fro — a m in Eq. (||) defines the (m- 
dependent) time when the measurement is started, and 
may not be easily observed in a transport experiment. 
To illustrate our point that I m (where m is an odd num- 
ber) disappears for certain orientations it is enough to 
plot the amplitude C m (V). 

We now investigate the ac components. For s-wave 
superconductors we have the ac Josephson effect where 
the Josephson frequency is ujj. In the case of <i-wave su- 
perconductors this may not be true, as discussed above. 
For the case d x 2_ y 2 jd xy we find G\ = and C 2 ^ 0, 
meaning that the Josephson frequency for this configura- 
tion is 2a; j. Decreasing a 2 from 7r/4 (the case d xy ), the 
magnitude C m of the odd ac components increases from 
zero (and therefore reestablishes the ordinary Josephson 
frequency ujj) until a 2 reaches 0, in which case the sit- 
uation more resembles the s-wave case. To understand 
this effect we study the case d x 2_ y 2 /d xy . The expres- 
sion for the current involves an integration over both 
positive and negative injection angles 9. We may ex- 
plicitly calculate the current for positive and negative 
angles separately, and then add them to get the final re- 
sult. When changing the angle 9 to —9 we must change 
the phase 4>o between the superconductors by 7r, because 
A 2 (-0) = -A 2 (0) - A 2 (9)e™ (meaning O -» <fo + n). 
The effect of this is seen in the wave function in Eq. (Q), 
where we get an extra phase —inn. This phase, separated 
out together with <^>o, will pass through the calculation 
and finally show up in the expression for the current: 

I(V, t) = 5^(I m (V, 9 > 0) + I m (V, 9 < 0))e imujt 

m 

= Y, I m(V,9> 0)(1 + e iBMr )e ima ' jt , (9) 

m 

where the amplitudes I m (V,9 > 0) [I m (V, 6 < 0)] now 
are defined with an integration over positive (negative) 
angles only. The factor (1 + e wmr ) is zero for odd m and 
two for even m. This implies a doubling of the Josephson 
frequency for the case d x 2_ y 2 / d xy . 

We present the ac components for different configura- 
tions of the superconductors and three different barrier 
strengths Z. The plots are shown in Fig. ||, and are only 
for V > since one can show that J(— V) = —I(V). For 
the d x 2_ y 2/d xy junction we only plot the first non-zero 
amplitude C 2 . For the d x 2_ y 2 /d a2 ^/4 junction we plot 
the components C\ and C 2 since they are both non-zero. 

The results show that the amplitude C% decreases 
smoothly as the angle a 2 is increased from to tt/4. 



There is subharmonic gap structure (SGS) for the case 
of non-zero barrier strength at the same, voltages as in the 
s-wave case [see Figs. |^(b) and_0(c)].ta Angular averag- 
ing partly washes out the SGS.EJThere is no SGS for t 
ballistic case Z — at zero temperature [see Fig. |^(a)]. 
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FIG. 2. Here we show the first two Fourier components as 
a function of voltage for three different barrier strengths Z. 
The orientation of the left superconductor is fixed at «i = 0, 
while we vary the orientation of the right superconductor from 
a?2 = to a2 — 7r/4. The first component Ci is identically 
zero for a.2 = 7r/4 as explained in the text. Zero temperature 
is assumed. 

Considering the second Fourier component C 2 , we first 
note that this component is much smaller than C\ sim- 
ply because Ci corresponds to a higher order term in the 
expansion of the current in Eq. (0). The curves depicted 
in Figs. |(e) and |(f) again show SGS \pq SGS in the 
ballistic case Z = shown in Fig. |(d)]B When CV2 IS 
close to 7r/4 and Z is large [see Fig. g(f)], MGS clearly 
influence the curves of C 2 , resulting in an increase in 
C 2 for small voltage. MGS appear as resonances in the 
scattering states of Eq. (|4|) when the order parameter ex- 
perienced by the quasiparticle changes sign after normal 
reflection at the barrier. Because of the large density 
of states at the gap, a quasiparticle trajectory passing 
through the MGS as the one depicted in Fig. g will give 
rise to a large contribution to the current. The reso- 
nant trajectory originates from the gap edge only when 
3eV — Ai = Aq cos 29 as seen in Fig. H, meaning that 
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the upper limit for this particular process is Ao/3. 

Another quasiparticle trajectory passing through MGS 
is the one travelling back and forth within the gaps twice 
instead of six times as shown in Fig. |j| This process will 
give rise to a smaller current contribution to C2 than the 
one shown in Fig. || for the following reason. Calculating 
C2 we study energies Em and E2 n +A- For the process 
corresponding to eV = A cos 29, E 2n +i will be above 
the gap region, and therefore this process will be sub- 
ject to two additional non-perfect (out of gap) Andreev 
reflections which suppress the contribution. 
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FIG. 3. The trajectory above gives rise to a large amplitude 
of the second component C2 at small voltages, see Fig. |i|(f ) . 
There is a resonance when the trajectory starts from the gap 
edge of the left superconductor, hits MGS, and finally ends 
at the gap again. 

The condition eV = A± = Aocos2# for the process 
shown in Fig. || does not depend on A2 . This means that 
SGS corresponding to this process is present for all angles 
a2 but the amplitude is smaller when a 2 7^ 7r /4 since the 
resonant MGS do not occur at as many injection angles 9 
in this case. For a 2 = the structure is absent since there 
is no MGS in this case. Larger transmission (smaller Z) 
broadens MGS, making the effects of MGS less drastic, 
as seen when comparing Fig. ||(e) with Fig. §(f). The 
first component C\ shows no increase in amplitude at 
small voltages. Actually, there is similar SGS but this is 
smaller in amplitude than the ordinary background and 
does not dominate in the same way as for C2. 

Considering the case of two d xy superconductors, we 
get a phase shift of it in both superconductors when we let 
9 — > — 9, leading to no overall phase difference. Therefore 
the first ac component C\ does not vanish in this case. 

In experiments on high-T c d x 2_ y 2 / d xy grain bound- 
ary junctionsjlj Shapiro steps have been found at volt- 
ages corresponding to a Josephson frequency of 2ujj. As 
pointed out in Ref. [l^ this effect may be explained by 
assuming that parallel junctions are formed at the grain 
boundary. According to our results it is possible that the 
doubled Josephson frequency is due to d-wave symmetry 
of the order parameter. 

Very recently, it has been shown that there-is a nonzero 
ground state current parallel to the junctionEj in the sta- 
tionary case (zero voltage) for a s/d xy junction. Inspired 
by these investigations, we have used our formalism to 
calculate the parallel current for the non-stationary case. 
For junctions with two d-wave superconductors we find 
both time-independent (m = 0) and time-dependent 
(m =/= 0) currents in the general case. For the ^-direction 
current the even components disappear in the d x 2_ y 2 j d xy 
geometry while the odd components are non-zero. This 



happens because we get a minus sign between the two 
terms in Eq. (||) instead of the plus sign, making the 
x-direction and the y-direction currents complementary. 
This resembles the stationary case where the parallel cur- 
rent appears when the perpendicular current is zeroEl 

Also for the N/d a junction we find a time indepen- 
dent (m = 0; no ac in the NS case) current density in 
the y-direction at x = if the up/down symmetry of the 
junction is broken by the order parameter, which is the 
case when the angle a has a value between and 7r/4. For 
a = and a = 7r/4 the current density in the y-directi©a. 
is zero. The N/d a case has recently been discussed.Kll 
One notes that even for an asymmetric s-wave order pa- 
rameter there are current contributions in the y-direction 
for the NS case. 

In summary, we have calculated the ac components for 
voltage biased ci-wave junctions. We find that the odd 
ac components of the current in the ^-direction vanish if 
the orientations of the d-wave superconductors is a,\ = 
and a 2 = t/4. For this case the Josephson frequency is 
doubled, which has also been found in experiments on 
grain boundary junctions.Ej The reason for this effect is 
that the contributions from injection angles +9 and — 9 
have exactly the same magnitudes, but opposite signs. 
The angle average will then cancel the odd components. 
We have also discussed the possibility of current flow in 
the direction parallel to the junction. 
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